Formation of heavy quasiparticle state in two-band Hubbard model 
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A realization of heavy fermion state is investigated on the basis of two-band Hubbard model. By 
means of the slave-boson mean-field approximation, it is shown that for the intermediate electron 
density, n e = 1.5, the inter-band Coulomb repulsion U strongly emphasizes initially small difference 
between bands, and easily stabilizes integral valence in the lower band. As a result, a strong 
renormalization takes place in the lower band and the mixing strength between two bands, ft gives 
rise to a sharp peak at the Fermi level in the quasiparticle density of states, as that obtained in the 
periodic Anderson model. In contrast to a simple insight that the Hund's-rule coupling J reduces 
the characteristic energy, it turns out to be almost irrelevant to the renormalization for J < U. The 
required conditions are suitable for LiV2C>4, the first observed heavy fermion compound in transition 
metal oxide. 

PACS numbers: 71.27.+a, 71.10.Fd 



susceptibility ||, 7 Li 
thermal expansion |)J, 
To| , resistivity |l|] and so 



I. INTRODUCTION 

The recent discovery of heavy fermion behavior in 
L1V2O4 [0 uncovers the latent possibilities to explore 
Kondo physics in the lattice d-electron systems, which 
is restricted so far to /-electron systems containing lan- 
thanide or actinide atoms fg||. The heavy fermion be- 
havior has been widely observed in various measure- 
ments, such as specific heat 
and 51 V NMR §§, fiSR @J 
quasielastic neutron scattering 

on [ fi1^]l3| . The low-energy physics is characterized by the 
large mass enhancement in the specific heat coefficient, 
7 ~ 0.21 J/V-mol K 2 with the Sommerfeld- Wilson ratio 
i?w ~ 1.71. The characteristic temperature of Kondo or 
spin fluctuation is estimated as T* ~ 30 K. With elevated 
temperature, the magnetic susceptibility approximately 
follows the Curie- Weiss law, C j (T — 9), where the Curie 
constant C is consistent with a V +i spin 5=1/2 with g 
factor 2.23, and the negative Weiss temperature (0 = —63 
K) is familiar to /-electron heavy fermions. 

The several band structure calculations have been 
made at present [O |l7| . They have revealed that the 
octahedral coordination of the oxygen ions around the 
V atom causes the large splitting of d states into t2 g 
and e g orbitals. The partially filled ti g bands can be 
described roughly by V-V hopping, and they are well 
separated by the filled 0-2p-bands and the empty e g 
bands. Eyert et al. suggest the specific heat enhance- 
ment comes from spin fluctuations with the magnetic 
order suppressed by the geometric frustration |14|| . In 
similar context, spin fluctuation nearby magnetically un- 
stable point in Lii-zZn^ V2O4 is discussed by Fujiwara 
et al. 

On the other hand, an attempt to map onto the con- 
ventional periodic Anderson model (PAM) is made by 



Anisimov et al. |||||jT^]. In the realistic treatment of 
the trigonal symmetry crystal field, triply degenerate t2 g 
orbitals split into the non-degenerate A\ g and doubly de- 
generate E g representations of the D$d group. Their as- 
sertion is that due to the Coulomb interaction among 
d electrons, one electron of the d 1 ' 5 configuration is lo- 
calized into the A\ g orbital and the rest partially fills a 
relatively broad conduction band made from E g orbitals. 

The idea seems to resolve the enormous differences 
from the isostructural LiTi204, which has 0.5 d electrons 
per Ti atom. It shows the relatively T-independent Pauli 
paramagnetism in susceptibility and the superconducting 
state below T c = 13.7 K, which is well described by the 
BCS theor y [BOpl] ]. The related discussion is also made 
by Varma [fgf - ^ 

Nevertheless, it is not easy to for d electron be local- 
ized because of the much larger spatial extent of d or- 
bitals than of / orbitals, unless the Ai g orbital is located 
much deeper than the E g band. In other words, the use 
of the PAM has no solid ground contrary to a naive ex- 
pectation. Moreover, the intra-band Coulomb repulsion 
does not work effectively to make heavy mass since each 
electron favors to place in different orbitals rather than 
in the same ones. 

In this paper, we clarify how a heavily renormalized 
quasiparticle ground state is realized in the two-band 
Hubbard model, leading to somewhat different physical 
picture from that argued in Ref. |15|. By using the slave- 
boson mean-field approximation J23| , p4[ , we discuss the 
importance of the inter-band Coulomb repulsion to make 
one of the d electrons be technically localized and provide 
a situation similar to that given by the PAM. Since the 
geometric frustration is expected to prevent a long-range 
order, we restrict our attention to a paramagnetic ground 
state. In the same stand, we argue about LiTi2C>4 as the 
case of smaller electron density, n e = 0.5. 
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II. MODEL AND FORMULATION 



where the renormalization factors qe and q are given by 



Let us start with two-band Hubbard model, which may 
capture the broad E g band (called as A) with the width 2 
eV and the narrow A\ g band (called as B) with the width 
1 eV. The latter is located at 0.1 eV lower than the for- 
mer due to the trigonal distortion. The d-d Coulomb 
interactions and the Hund's-rule coupling are estimated 
as 3 eV and 1 eV, respectively |l5j . The mixing strength 
between bands must be much smaller than the width of 
both bands. The Hamiltonian is given by 



zj, q = VqlqE, 
1 



kia 



Hint = Utnu^nui + U ^Aa^Aa^Ba^Bfi 



J 



^2 4A a c *Af3c\ Bl c lB5 a al3 ■ a l5 . (1) 

The first term denotes the kinetic energy of conduction 
electrons for bands £, in which Ee = ±A/2 for A and B 
bands, A being the trigonal splitting, and fi the chemi- 
cal potential. The second is a mixing strength between 
bands and its fc-dependence is neglected for simplicity. 
The intra- and inter-band Coulomb interactions Ue and 
U as well as the Hund's-rule coupling J at the site i are 
considered in Hi nt , where the intra-band Coulomb inter- 
actions are set by XJ a — Ub = oo for simplicity. 

To solve this Hamiltonian, we introduce slave-boson 
fields for d°-d 2 states at each site. We associate a bo- 
son a for d° state, pa a for d 1 states, and diss* f° r d 2 
states labeled by their spin states, (S,S Z ), respectively. 
For the uniform solution of the mean-field approximation, 
we replace these bosons by site- independent c numbers. 
Assuming the paramagnetic ground state, five bosons, 
e = e», pa = PiAa, Pb = PiBa, d = d 0Q and di = dis z , 
are involved in calculation. 

The completeness relation for dP-d 2 states is given by 



1 = 0, 



I = e 2 



d 2 



3d 2 . 



(2) 



Since the probabilities for the singlet and the triplet 
states are given by dp an d 3d 2 , respectively, the two-body 
interactions are rewritten in terms of bosons as 



H^/N = U(d 2 + 3d 2 ) - -J (d? - dl) 



(3) 



where N is the number of sites. The each species of elec- 
trons must satisfy the constraint at each site, 



Qz = 0, Q e 



i 



(d 2 +3d 2 ) 



(4) 



In this slave-boson scheme [ p3| , the hopping term is 
suppressed as 



qc\ Aa c kB< y, (5) 



zt = e p e + - (d + 3di) p- t . 



£ - 1/2 ^(l-Q,)- 1/2 , (6) 
(7) 



Note that the Gutzwiller correction Q e ' (1 — Q 
necessary to reproduce non- interacting limit [^3 25 1. 

Finally, we obtain the mean-field free energy per site 
with two Lagrange multipliers A and Xe, 



-1/2 



+ ffint, F mf /n = _JLj2J2 ln ( 1 + e-^*"-")) + H^/N 



k m 



+X(I- 1) -2j2^eQe, 



(8) 



where the bonding and the anti-bonding (m = =p) quasi- 
particle bands are given by 



Ekm — 



1 



£,ki = qieke + Ee + Xt, v = qv. 



(9) 



The width of the band I and the mixing strength are 
renormalized by factors qe and q, respectively. The posi- 
tion of the bands are moved up by an amount A^. 

Minimizing the free energy with respect to five bosons 
and two Lagrange multipliers, the set of self-consistent 
equations are obtained as follows: 



+ 2]>> 2 



3d\ 



1 = 0, 



n i -2\pj+-(d 2 + 3di)} = 0, 

vf^- + Xe — 0, 

oe 



EQ dqi 
2 de 

i 

Q dqe 
2 dpe 

Q dqe 
2 dd 

Q dqe 



E 

E 

i 

E 



_ dq 

vr- 2(Xe> - X)pe = 0, 

dpi' 



vr 



■ vr- 



dq 
dd 

dq 



2 ddi ddi 

' E fit = °' 



(T - A ^ + A ) rf o = 

e 

3(Ti - J2 ^ + A ) rf i 



0. 



(10) 

(11) 
(12) 

(13) 

(14) 

(15) 
(16) 



where the energies of the singlet and the triplet states 
are defined by T = U + 3 J/4 and T x = U - J/4. The 
last equation is responsible for determining the chemical 
potential for given electron density n e . Hereafter we re- 
strict ourselves to the case at zero temperature /3 _1 = 0. 
At zero temperature the averages of electron densities, 
mixing amplitude and kinetic energies are given by 
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Nnt = {At^ki.) = £ pftOQji - E km \ (17) 

ka kma 

Nf = E Ua^b.) = E cr^ - 4m), (is) 



t e - 1 l-n B 



(22) 



Nee = E e " ( c L<r c w<»-) = E £ ^/ 9 S^(m ~ Sfem), (19) 



with 



Pfcf — 2 



1 + m 



(aA-as) 2 +45 2 



(6a-6s) 2 +4w 2 



(20) 
(21) 



For simplicity, we use a rectangular density of states 
(DOS) with a linear dispersion relation, i.e., eke = Wex/2 
for \x\ < 1. Then, The fc-summation in the aver- 
ages can be carried out analytically with the integration, 



III. RESULTS AND DISCUSSIONS 



The set of self-consistent equations are solved numer- 
ically. In the following we use parameters, Wa = 2, 
Wb = 1, and v = 0.2 (eV). Figure [I] shows the quasi- 
particle DOS with and without interactions, Ue, U and 
J, for A = 0.2 (eV) and n e — 1.5 corresponding to 
L1V2O4. The bandwidth is renormalized slightly by the 
intra-band Coulomb repulsion Ua = Ub = 00. However, 
the renormalization amplitude of the narrower B band, 
qB, remains at the order of 10 _1 without the inter-band 
repulsion, since each electron favors to place in different 
orbitals rather than in the same ones. On the contrary, 
with the inter-band interactions U — 3 and J = 1 (eV), 
a strong renormalization takes place and it gives rise to a 
sharp peak at the Fermi level in the quasiparticlc DOS. 
(Its width is about 40 K.) Note that both the upper and 
the lower Hubbard bands cannot be argued in the mean- 
field approximation. The inset in Fig. fu shows the quasi- 
particle DOS for the case of n e = 0.5 corresponding to 
LiTi204. As it is expected, the renormalization is very 
weak, qs ~ 0.7. 

To elucidate why the inter-band interaction assists the 
strong renormalization, we discuss the limiting cases for 
n e = 1.5. In the absence of Hi nt in eq. (0), a rather 
large trigonal splitting, i.e., A ~ is required to 

stabilize the integral valence in the B band, ub- While, 
in the case of the strong repulsion, Ua — Ub — 00, and 
U/We 1, the inter-band Coulomb repulsion consid- 
erably enhances the difference of electron densities be- 
tween two bands, An = Tib — nA, because of the relation 
UfiAUB = U[n 2 e — (An) 2 ]/4. In this case with do ~ d\ 
for J/U <C 1, the renormalization factor vanishes as 



Note that in the PAM with U = 00 the hybridization be- 
tween / and conduction electrons is suppressed as V 2 — > 
q f V 2 , where q f = (1 - n/)/(l - n f /2) |i||§j2|. We 
emphasize here that although the mechanism of strong 
renormalization appears similar to that in the PAM, it is 
totally a new mechanism that the Kondo limit is dynam- 
ically provided by the inter-band Coulomb repulsion. It 
will be shown below that the integral Tib can be stabilized 
even for rather small A, which is based on the detailed 
balance between the kinetic energy and the inter-band 
Coulomb interaction. 

On the other hand, in the case of J ^ (U, We), the 
amplitude of triplet state, d\, becomes as large as possi- 
ble. Its maximum is bounded by 3d 2 < min(n,4,nB). 
Thus, in order to take the largest value, An is sup- 
pressed. In the limit of large J, the probabilities for 
dP-d 2 states, (e 2 ,p 2 ,3d 2 ) approach (1 — n e /2, 0, n e /2), re- 
spectively. Namely, the system undergoes a dimerization 
with a charge order, and the renormalization factor qg 
vanishes since pi — > in eq. (Q). 

Figure || shows Tib as a function of the trigonal split- 
ting A. The intra-band Coulomb interaction Ua and Ub 
(circle) somewhat enhances To stabilize the integral 
valence, however, A is required as large as that for the 
case of free electrons (square), i.e., A ~ 0.5 eV. On the 
other hand, the inter-band interaction (triangle) works 
effectively to stabilize integral valence even for small A. 
Note that ub is almost unity for A ~ 0.1 eV. 

The A dependence of the renormalization factors are 
shown in Fig. ||. As expected from the above discussion, 
in the presence of the inter-band interactions (square), B 
band is highly renormalized owing to hs — > 1, while the 
intra-band interactions are almost irrelevant up to rela- 
tively large A (circle). It is noted that the upper band 
is almost unrenormalized, while the band mixing is also 
renormalized considerably, i.e., q = ^JqAqB- 

In order to elucidate how the inter-band interactions 
reduce renormalization factors, we extract the interaction 
dependences of the renormalization factors for A = 0.2 in 
Fig. |J. It is shown that the inter-band Coulomb interac- 
tion (square) effectively reduces renormalization factors. 
On the other hand, the Hund's-rule coupling turns out 
to be almost irrelevant for J < U (circle, triangle). This 
is in contrast to a simple insight that a strong cancel- 
lation between the Hund's-rule coupling and the Kondo 
exchange coupling considerably reduces the characteris- 
tic energy as discussed in Refs. J2^,^§]. Discussions are 
based on the impurity model and hence there is no con- 
straint for the electron density such as n e = ua + nB- 
Since a change of Hund's-rule coupling requires another 
change of parameters to restore a given electron density, 
all parameters must be treated in a self-consistent fash- 
ion. This might remove the discrepancy from our re- 
sult that the Hund's-rule coupling is almost irrelevant in 
renormalization. 



3 



IV. SUMMARY 

We have shown that the inter-band Coulomb repulsion 
plays a significant role to reduce the renormalization fac- 
tor strongly, while the Hund's-rule coupling is almost ir- 
relevant in renormalization for J < U . Even though both 
bands are rather broad and a splitting between bands is 
very small, the resultant quasiparticle can have a heavy 
mass enhanced by about 10 2 times since the Kondo limit 
is dynamically provided by the inter-band Coulomb re- 
pulsion. Although the values obtained by slave-boson 
approach may be changed quantitatively by more elabo- 
rate one, the situation is highly plausible to account for 
the heavy- fermion behavior in LiV204 and the enormous 
differences from LiTi204. 

Since the heavily renormalized quasiparticle has been 
stabilized dynamically by the inter-band Coulomb repul- 
sion, it should couple strongly with orbital fluctuations 
at higher temperature. The large contribution to the spe- 
cific heat observed above T* j|] is presumably related to 
the orbital fluctuations. 

At low temperature d electron systems generally ex- 
hibit a long-range order. If a paramagnetic state survives 
due to some reasons such as a geometric frustration, one 
would expect heavy-fermion behavior in numerous d elec- 
tron systems. The resultant quasiparticle holds the pos- 
sibility of showing fascinating phenomena such as a novel 
superconductivity mediated by orbital fluctuations. 
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FIG. 3. The A dependence of renormalization factors 
for Wa = 2,Wb = 1,v = 0.2 (eV) and n e = 1.5. The B 
band is highly renormalized by the inter-band interactions. 



FIG. 1. The quasiparticle density of states with and 
without interactions for Wa = 2, Wb = 1, v = 0.2, 
A = 0.2 (eV) and n e = 1.5 corresponding to L1V2O4. 
The inset represents the case of L1T12O4 as n e = 0.5. 
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FIG. 4. The inter-band interaction dependences in the 
renormalization factors for Wa = 2, Wb = 1, v = 0.2, 
A = 0.2 (eV) and n e = 1.5. 



FIG. 2. The A dependence of electron density in B band 
for Wa = 2, W B = 1, v = 0.2 (eV) and n e = 1.5. The 
comparison among three cases shows that the inter-band 
interactions effectively works to stabilize integral valence 
even for rather small A. 
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